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Abstract 

This  report  summarises  results  of  coalescence  efficiency  of 
colliding  drops  which  were  obtained  from  laboratory  experiments  and 
a theoretical  model.  Both  theory  and  experiment  suggest  that  the 
coalescence  efficiency  is  a function  of  both  drop  radii  and  drop  size 
ratio.  Even  though  no  clear  dependence  of  coalescence  efficiency  on 
impact  speed  was  observed  experimentally,  the  theory  does  predict  some 
dependence  on  velocity.  Generally  it  is  observed  that  the  coalescence 
efficiency  decreases  with  increase  in  the  size  of  the  large  drop,  Rj^, 
and  the  drop  size  ratio,  p.  However,  for  small  a decrease  in  the 
efficiency  was  also  shown. 

Both  theory  and  experiment  demonstrate  the  strong  dependence  of 
coalescence  on  impact  angle.  The  theory  is  able  to  predict  the  efficiency 
of  coalescence  for  different  Impact  speeds  and  not  only  the  relative 
terminal  velocities  of  the  drops.  It  is  also  capable  of  calculating 
the  coalescence  efficiency  of  all  impacts  on  stationary  targets  and  the 
coalescence  efficiencies  for  moving  drops  when  the  p ratio  is  small. 

From  the  experiments  two  empirical  equations  were  obtained  which  relate 
the  coalescence  efficiencies  and  partial  coalescence  efficiencies  with 
p and  R^.  These  could  easily  be  used  in  numerical  models  of  cloud  and 
rain  development. 


CHAPTER  1:  INTRODUCTION 


The  growth  of  warm  clouds  is  understood  to  take  place  by  three  dis- 
tinct processes:  a)  the  nucleation  process,  whereby  droplet  embryos  are 

formed  on  aerosol  particles  (called  cloud  condensation  nuclei),  b)  the 
growth  process  by  condensation  which  is  the  process  of  diffusion  of  vapor 
from  the  environment  to  the  drops,  and  c)  the  growth  by  the  collision- 
coalescence  process  in  which  drops  grow  by  capturing  smaller  droplets 
during  their  fall  or  rise  in  the  cloud.  It  was  believed  that  the  growth 
by  condensation  dominates  the  growth  of  droplets  up  to  a radius  of  about 
20y,  while  the  collision-coalescence  process  dominates  the  growth  of  larger 
drops.  Fog,  on  the  other  hand,  grows  initially  by  the  nucleation  process 
and  continues  to  grow  by  condensation.  Since  very  little  convection 
exists,  the  fog  droplets  tend  to  remain  small  and  their  growth  by  collision- 
coalescence  is  negligible.  However,  some  of  the  artificial  fog  dispersal 
attempts  have  used  different  techniques  such  as  the  creation  of  convection, 
introduction  of  steam,  spread  of  hygroscopic  particles  and  the  spraying 
of  highly  charged  droplets.  The  dispersal  of  the  fog  is  then  accomplished 
by  either  evaporation  of  the  droplets  by  mixing  them  in  a subsaturated 
environment  or  by  forcing  the  droplets  to  grow  to  such  sizes  that  they 
precipitate  to  the  ground.  The  latter  growth  is  achieved  either  by  con- 
densation of  excess  water  vapor  onto  the  existing  droplets  or  by  coalescence 
of  the  fog  droplets  with  the  artificially  produced  ones. 

A great  deal  of  theoretical  work  has  been  done  on  the  collision 
efficiencies  of  water  drops  in  the  atmosphere  (e.g.  Hocking,  1959,  Shafrir 
and  Neiburger,  1963,  and  Hocking  and  Jonas,  1970).  All  the  results  indicate 
that  the  larger  the  drop  gets,  the  higher  is  its  collision  efficiency  for 
the  same  p = r/R  ratio.  As  p increases  (for  the  same  large  drop  size)  from 
0 to  0.5  the  efficiency  increases.  For  larger  p the  efficiency  drops  off 
only  to  rise  again  at  p closer  to  unity. 

The  number  of  experiments  testing  these  theories  have  been  limited 
in  number  and,  as  a result  of  experimental  difficulties,  also  limited  to 
very  low  p ratios  ( <0.3)  and  to  the  very  high  ones  ('V'1.0)  (e.g.  Picknett, 
1960,  Woods  and  Mason,  1954,  Beard  and  Pruppacher,  1971).  In  these  experi- 
ments the  collision  efficiency  could  not  be  measured  independently  and 
could  only  be  deduced  from  measurements  of  collection  efficiency  by  assuming 
a unity  coalescence  efficiency.  Since  their  results  agreed  fairly  well  with 
the  theoretical  predictions  of  the  collision  efficiencies,  they  concluded 
that,  in  the  range  they  have  tested,  the  coalescence  efficiency  is  indeed 
unity. 

Most  numerical  models  for  the  growth  of  cloud  droplets  or  the  dis- 
persal of  fog  by  introducing  larger  drops  have,  for  lack  of  better  values. 


assumed  that  the  coalescence,  is  unity.  Recently  Neiburger  et  al.  (1972) 
and  Levin  et  al.  (1973)  have  reported  results  of  son*  laboratory  experi- 
ments in  which  drops  of  various  sizes  were  suspended  in  moving  clouds  of 
smaller  droplets.  The  liquid  water  content,  the  s<ze  distribution  of 
the  smaller  droplets,  the  temperature  and  saturation  ratio  were  carefully 
controlled  and  monitored.  Their  results  indicated  that  drops  smaller 
than  40y  in  radius  did  not  grow  at  all  by  the  collision-coalescence  pro- 
cess after  being  suspended  in  a cloud  of  droplets  of  mean  radius  of  12y. 
Larger  drops  ( > 55y)  grew  by  collecting  smaller  droplets  (p''^0.2)  but  at 
a much  slower  rate  than  predicted.  When  electric  charges  were  applied 
to  the  drops,  their  rate  of  growth  Increased. 

Their  conclusions  were  that;  a)  the  coalescence  efficiency  cannot 
be  considered  unity  for  all  p ratios;  b)  that  as  the  p ratio  increases 
(at  least  to  p = 0.35)  the  coalescence  efficiency  decreases.  These 
results  have  grave  consequences  on  the  theory  of  growth  of  warm  clouds 
or  the  effectiveness  of  the  collision-coalescence  mechanism  in  the  arti- 
ficial dispersion  of  warm  fogs,  since  if  the  coalescence  efficiency  is 
lower  than  unity,  it  is  difficult  to  understand  how  these  clouds,  for 
exauiple,  can  grow  at  all.  In  fact  some  doubt  has  already  been  cast  on 
the  effectiveness  of  the  collection  mechanism  process  to  account  for 
the  rapid  growth  of  warm  clouds  even  with  coalescence  efficiency  of 
unity  (see  Bartlett,  1970).  Since  the  results  of  Levin  et  al.  (1973) 
indicate  even  lower  coalescence  efficiencies,  it  is  even  more  difficult 
to  explain  the  growth  of  these  clouds  by  this  mechanism. 

However,  their  conclusions  were  also  based  on  the  observations  of 
the  collection  and  not  on  an  independent  investigation  of  the  coalescence. 
They  assumed  that  the  theoretically  computed  collision  efficiencies  that 
had  been  previously  proposed  were  correct,  and  from  the  ratio  of  their 
experimentally  obtained  collection  efficiencies  and  the  theoretical  colli- 
sion efficiencies,  they  computed  the  coalescence  efficiencies 
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Only  a few  reported  experimental  results  are  available  in  which  the 
coalescence  efficiency  was  independently  investigated  (e.g.  Charles  and 
Mason,  1960).  Recently  Whelpdale  and  List  (1971)  have  reported  results 
of  an  experiment  in  which  coalescence  on  drops  with  p ratios  of  the  order 
of  0.1  were  investigated.  Although  the  sizes  they  used  and  the  p ratios 
apply  to  rain  drops,  their  results  indicate  somewhat  the  same  trends  in 
the  coalescence  efficiencies  as  was  found  by  Levin  et  al.  (1973).  They 
found  that  their  results  could  be  described  by  the  equation 
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However,  since  their  measurements  were  carried  out  in  a narrow  range  of 
p ratios  and  only  with  a few  drop  sizes  the  above  equation  could  be  con- 
sidered as  an  extrapolation.  Braize r-Smith  et  al  (1972)  investigated  the 
separation  conditions  for  interacting  drops  and  the  number  of  fragments 
so  produced.  However,  their  results  were  limited  only  to  p ratios  between 
0.4-1  and  as  a result  of  their  experimental  set  up,  did  not  eliminate 
the  flow  around  both  interacting  drops. 

Our  aims  in  this  work  were:  (1)  To  experimentally  extend  the  range 

of  measurements  carried  out  by  Whelpdale  and  List  (1971)  to  include  a 
variety  of  drop  sizes  and  drop  ratios.  (2)  To  obtain  an  empirical  equation 
to  describe  the  coalescence  and  partial  coalescence  efficiencies  for  use 
in  numerical  models  of  cloud  growth,  electrical  development  and  fog  dis- 
persal. (3)  To  theoretically  investigate  some  of  the  main  forces  hindering 
the  coalescence  process  and  mathematically  describe  the  motion  of  approaching 
drops.  The  main  point  in  this  theoretical  investigation  is  the  effect  of 
the  drop  deformation  on  the  coalescence  process. 


CHAPTER  2:  THE  COALESCENCE  EXPERIMENT 

2-L  THE  EXPERIMENT 

In  order  to  separate  the  problem  of  coalescence  from  that  of  the  colli- 
sion we  used  a similar  approach  to  that  of  Whelpdale  and  List  (1971).  In 
this  experiment  large  drops  were  suspended  from  a small  capillary  tube.  A 
small  droplet  generated  by  a drop  generator,  impacted  the  larger  drop  on 
its  lower  hemisphere.  Viewing  the  interactions  from  two  perpedicular 
directions  permitted  the  identification  ot  theangle  of  impact.  The  inter- 
actions were  simultaneously  photographed  by  a 35mm  camera  and  by  a television 
camera.  This  method  of  viewing  was  thought  to  be  an  in^rovement  over  the 
method  used  by  Whelpdale  and  List  who  observed  the  interaction  from  a single 
direction  and,  hence,  could  have  had  difficulties  in  determining  the  exact 
inpact  angle.  As  will  be  shown  later  the  impact  angle  is  an  important  para- 
meter for  the  coalescence  and  uncertainty  in  its  determination  could  lead  to 
great  inaccuracies.  In  order  to  separate  the  collision  from  the  coalescence 
problem  one  has  to  eliminate  the  flow  around  the  interacting  drops.  This 
is  an  impossible  task.  However,  the  method  used  by  Whelpdale  and  List  and 
also  used  here  seems  to  be  a close  approximation  to  it.  When  a small  drop 
approaches  a larger  falling  drop  it  enters  the  region  of  the  boundary  layer 
of  the  larger  drop.  Within  this  region  the  flow  decreases  as  one  approaches 
the  drop  surface.  For  large  drops  the  flow  does  not  diminish  at  the  surface 
due  to  internal  circulation  in  them.  However,  as  the  drops  approach  each 
other  the  flow  between  them  is  probably  controlled  by  the  drainage  of  the 
intervening  air  layer.  This  drainage  depends  on  the  speed  of  the  approaching 
drops  and  their  deformation.  Therefore  it  seems  that  the  experiment  of  Whelpdale 


and  List  and  the  one  proposed  here  are  applicable  to  interaction  of  drops  in  the 
atmosphere,  at  least  for  small  p ratios.  At  larger  p ratios  and  larger 
drops  it  is  possible  that  the  leading  surface  of  the  smaller  droplet  enters 
the  boundary  layer  while  the  other  side  is  still  affected  by  the  air  flow 
further  out.  In  such  a situation  the  flow  will  tend  to  move  the  droplet 
sideways  and  reduce  the  probability  of  coalescence.  Hence  for  such  cases 
the  results  that  we  obtain  here  can  be  considered  as  an  upper  limit  to  the 
coalescence  efficiency. 


2-2.  THE  APPARATUS 

Fig.  2-1  is  a schematic  diagram  of  the  experimental  set  up.  The  size 
of  the  larger  collector  drop  was  controlled  by  a micrometer  syringe  and 
could  be  varied  from  200  microns  to  over  2.5mm  in  diameter.  Its  position 
in  space  could  be  changed  by  a micromanipulator.  The  smaller  droplet  was 
generated  by  a drop  generator  similar  but  more  powerful  than  that  of  Abbott 
and  Cannon  (1972).  The  generator  consists  of  an  acoustic  suspension  speaker 
with  a metal  rod  connected  to  the  center  of  its  cone.  The  other  side  of 
the  rod  is  bent  and  immersed  in  a water  bath.  An  electrical  pulse  of  known 
amplitude  and  shape  pushes  the  speaker's  cone  up  and  forces  the  rod  out  of 
the  water  and  produces  a droplet.  The  droplet  size  is  controlled  by  the 
size  of  the  rod,  the  pulse  shape  and  its  amplitude.  It  was  constructed 
in  such  a way  as  to  be  able  to  eject  vertically  upward  single  small  droplets. 
The  generator  was  directed  with  a slight  angle  to  the  vertical  and  therefore 
ejected  droplets  in  a parabolic  trajectory.  This  trajectory  prevented  the 
droplets  from  falling  back  into  the  water  reservoir  and  disturbing  the  next 
generated  droplet.  By  positioning  the  large  drop  at  different  heights  along 
the  upward  trajectory  of  the  droplet  different  impact  velocities  were  obtained. 
By  horizontally  varying  the  position  of  the  collector  drop  the  impact  angles 
could  be  varied.  Each  interaction  was  simultaneously  photographed  on  a 35mm 
film  and  recorded  on  a video  tape  by  a television  camera  which  was  mounted 
at  a right  angle  to  the  viewing  direction  of  the  still  camera.  The  inter- 
ation  zone  was  illuminated  by  a continuous  flood  light  (lOOOw)  with  a 
water  filter  in  front,  and  by  a stroboscope  flashing  at  exactly  1000  cycles 
per  sec.  The  position  of  the  flood  light  and  the  stroboscope  were  at  160° 
and  180°  from  the  direction  of  observation  of  the  35mm  camera. 

Once  the  position  of  the  trajectory  was  determined  the  electronic 
circuit  was  triggered.  The  drop  generator  was  the  first  to  be  triggered 
and  was  followed  after  pre-set  delays  by  the  opening  of  the  camera's 
shutter  and  ten  flashes  of  the  strobe.  Then  the  camera  shutter  was  closed 
automatically  and  the  film  advanced.  During  this  whole  period  the  television 
camera  continuously  recorded  the  entire  event. 
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By  adjusting  the  position  of  the  inpact  place  so  that  collisions 
would  occur  at  precisely  the  vertical  great  circle  of  the  collector  drop 
and  in  the  focal  plane  of  the  camera,  the  exact  angle  and  drop  size  could 
be  determined  from  the  photographs.  The  impact  velocity  was  determined 
from  the  distance  between  successive  light  spots  on  the  film.  The  exact 
trajectory  was  determined  from  the  streak  produced  by  the  droplet  due  to 
the  continuous  illumination  by  the  flood  light. 

Both  the  collector  drop  holder  and  the  droplet  generator  were 
grounded.  This  eliminated  electrical  charges  from  the  collector  drop  but 
not  from  the  droplets.  However,  it  was  estimated  that  the  maximum  charge 
they  carried  did  not  exceed  10“  esu  and,  hence,  was  believed  to  have 
little  effect  on  the  coalescence  (eg.  Neiburger  et  al,  1974). 


2-3.  RESULTS 

About  8000  photographs  were  taken  out  of  which  over  3000  were  found 
adequate  for  analysis.  Collisions  studied  could  be  classified  into  three 
types:  (1)  those  resulting  in  full  coalescence,  (2)  those  resulting  in 

bouncing  and  (3)  those  resulting  in  partial  coalescence.  The  difference 
between  the  first  and  the  last  two  types  of  interactions  could  easily  be 
identified  from  the  photographs.  However,  the  differences  between  the 
last  two  types  was  more  difficult  to  resolve.  Partial  coalescence  events 
where  identified  as  such  when  the  size  of  the  outgoing  droplet  was  smaller 
than  its  original  size  before  impact.  However,  since  differences  of  only 
about  20  microns  in  diameter  could  be  resolved,  many  Interaction  which 
actually  resulted  in  partial  coalescence  may  have  been  classified  as  bounce. 
This  may  mean  that  the  efficiency  of  partial  coalescence  we  obtained 
actually  represents  a lower  bound. 

Fig.  2-2a,b  and  c represent  interactions  which  resulted  in  coalescence, 
bounce,  and  partial  coalescence  respectively  (details  are  given  in  the 
figure  captions). 

In  order  to  resolve  the  possible  effect  of  the  impact  angle  and  impact 
velocity  on  the  coalescence  and  bouncing  events  all  the  results  were  gathered 
into  two  graphs  Fig.  2-3a  and  2-3b.  In  here  no  reference  is  made  to  the 
drop  size  or  drop  size  ratio.  It  is  observed  that  up  to  a velocity  of  2.4 
m/sec  there  is  no  clear  dependence  on  the  Impact  velocity,  however,  a depen- 
dence on  the  impact  angle  is  observed.  Coalescence  occurs  at  angles  close 
to  head-on  collision  and  decreases  as  the  angle  increases.  The  reverse  is 
true  for  bouncing.  However,  an  in?)ortant  conclusion  can  already  be  drawn: 
bouncing  is  observed  at  very  low  Impact  angles  in  contrast  to  observations 
by  Whelpdale  and  List  (1971). 
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Fig.  2-2a.  Coalescence  event  of  two 
drops  2200  and  110  microns  diameter 
impacting  at  a velocity  of  Im/sec 
and  angle  of  38°. 


Fig.  2-2b.  A bouncing  of  two  drops 
2180  and  130  microns  diameter  impacting 
at  a velocity  of  l.Om/sec  and  an  angle 
of  8°. 


Fig.  2-2c.  A partial  coalescence 
event  when  drops  of  2160  and  120 
microns  diameter  impact  at  a velocity 
of  0.5m/sec  and  an  angle  of  29°.  Afte 
separation  the  smaller  droplet  has  a 
diameter  of  80  microns  and  a velocity 
of  0.34m/sec. 
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Flg.2-3b.  Ex;^orlr',ental.  observation  of  bcnnclnp,  (•)  and  partial  Coaloacfnce  ("*) 
function  of  Velocity  imd  a-;gle  of  cc*nt;)ct. 
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It  is  Importaat  to  check  the  dependence  of  the  coalescence  frequency 
on  the  Ijnpact  angle.  Fig.  2-4  represents  the  coalescence  frequency  as  a 
function  of  inpact  angle  for  four  equal  angle  intervals  between  0°-90°. 

It  is  seen  that  no  clear  consistent  dependence  of  the  coalescence  frequency 
is  observed  as  a function  of  the  impact  velocity. 

Fig.  2-5  represents  the  coalescence  frequency  as  a function  of  inpact 
velocity  for  different  p ratios.  Again  no  clear  consistent  dependence  on 
inpact  speed  is  observed  although  some  decrease  in  coalescence  with  an  in- 
crease in  impact  speed  can  be  seen. 

This  lack  of  consistent  trend  in  the  dependence  of  the  coalescence 
frequency  on  the  impact  speeds  (of  up  to  2.4  m/sec)  led  us  to  the  idea  of 
gathering  the  data  without  regard  to  inpact  speed  and  instead  try  to  find 
some  dependence  on  p ratio,  impact  angle  and  the  size  of  the  larger  drop. 

Fig.  2-6  presents  the  coalescence  frequency  as  a function  of  impact 
angle  for  different  p ratios.  In  this  graph  no  separation  of  the  size  of 
the  large  drops  was  made.  The  conclusion  that  can  be  drawn  from  this  graph 
is  that  the  coalescence  decreases  rapidly  with  increasing  impact  angle  and 
with  increase  in  p ratio.  This  conclusion  is  in  agreement  with  the  results 
of  Whelpdale  and  List  (1971). 

The  last  possible  dependence  that  should  be  resolved  is  the  dependence 
on  the  size  of  the  larger  collector  drop. 

Figs.  2-7a  to  2-7f  represent  the  dependence  of  the  coalescence  frequency 
on  the  impact  angle  for  a variety  of  drop  sizes  and  p ratios.  As  expected, 
the  dependence  of  the  coalescence  on  the  Impact  angle  and  p is  observed  but 
in  addition  some  dependence  on  the  drop  diameter,  a,  also  emerges. 


2-4.  COALESCENCE  EFFICIENCY 

The  results  so  far  were  expressed  in  terms  of  coalescence  frequency. 
This  means:  the  fraction  of  events  resulting  in  coalescence  out  of  the 
total ^number  of  Interactions  taking  place  within  the  angle  interval  (eg. 
0-22.5).  Since  in  confutation  of  cloud  growth  it  is  easier  to  deal  with 
probabilities  or  efficiencies  of  coalescence  rather  than  frequencies  the 
above  results  have  to  be  converted.  For  this  purpose  we  adapted  the  method 
used  by  Whelpdale  and  List  (1971),  The  contribution  to  the  coalescence 
efficiency  of  interacting  drops  in  each  angle  interval  is  evaluated  by 
the  product  of  the  coalescence  frequency  in  this  interval  and  the  projected 
area  of  the  interval  onto  the  horizontal. 
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Fig.  2-4.  X of  coalescence  frequency  for  different  velocities  (n/sec.) 
of  small  droplets  for  Various  Contact  angles. 
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Fig.  2-r7c.  The  dependence  of  the  coalescence  frequency  on  the 
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The  projected  area  can  be  expressed  as: 


A - TTR^[8in^9„  - sln^ej 
B A 

where  and  0„  are  shown  in  Fig.  2-8  and  R is  the  radius  of  the  collector 
, A B 

drop. 


Figs.  2-9a-f  represent  the  histograms  which  are  the  result  of  the 
product  of  the  area  and  the  corresponding  coalescence  frequency  in  each 
angle  interval.  The  areas  under  the  histograms  represent  the  coalescence 
efficiencies.  Fig.  2-10  is  a representation  of  the  coalescence  efficiency 
as  a function  of  the  drop  size  ratio  (p)  for  different  large  drop  sizes. 

In  the  figure  a-  “ 200-300y  and  so  on.  It  is  clearly  observed  that  the 
coalescence  efficiency  decreases  as  the  p ratio  increases  as  was  observed 
by  Whelpdale  and  List  (1971)  and  Levin  et  al  (1973)  but  contrary  to  their 
results  a strong  dependence  on  the  size  of  the  large  drop  is  also  observed. 
Since  the  coalescence  is  found  to  be  a function  of  two  parameters  we  obtained 
a best  fit  surface  to  the  results.  This  surface  is  shown  in  Fig.  2-11  and 
can  be  represented  by  the  following  empirical  equation: 

E - 0.65  - 2.08x  10-3  RR  + 6.56  xlO-"*  p^R  + 6.07 
X 10-7  pr2  + i.  17  X 10-2  p2  - 9.3  x10-®r2 


where  R is  in  microns.  In  the  figure  the  dots  represent  the  experimental 
values.  The  mean  radii  were  chosen  so  that  for  exaple  a»  is  represented  by 
a point  at  r - a2/2  = 125y  and  similarly  for  the  other  a^s. 


Since  many  of  the  interactions  resulted  in  partial  coalescence  (events 
in  which  the  incoming  droplet  lost  part  of  its  mass)  it  is  valuable  to  analyze 
these  data  in  a similar  way  to  that  discussed  above.  Unfortunately,  out  of 
the  many  photographs  analyzed  not  enough  data  were  obtained  to  be  able  and 
separate  them  into  categories  based  on  both  large  drop  size  and  drop  size 
ratio.  Instead,  the  data  were  grouped  in  classes  according  to  p ratios  only. 
Fig.  2-12  represents  the  partial  coalescence  efficiency  in  each  angle  interval 
for  different  p ratios.  The  most  striking  feature  of  this  graph  is  the  fact 
that  in  most,  cases  the  highest  efficiency  for  partial  coalescence  is  observed 
between  45°  and  67.5°  even  though  a small  probability  for  such  occurances  is 
found  at  other  impact  angles, 
expressed  empirically  by: 

E 

partial  coalesce 
and  Fig.  2-13  represents  it  in 


The  equation  of  best  fit  to  this  data  can  be 

= 0.22  - 0.23p  - 0.025p2  (2-2) 

a graphical  form. 
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FI3.  2-9a.  Coalescence  efficiencies  a.s  a function  of  lipact  a.ngle 
for  d'ffereut  p ratios  end  sizes  of  large  drops. 
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Fig.  2-9c^  Coalescence  efficiencies  as  a function  of  Impact  angle 
for  different  p ratios  and  $lzes  of  large  drops. 
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2-5,  THE  EFFECTS  ON  THE  DEVELOPMENT  OF  THE  RAIN  DROP  SPECTRA 

The  above  results  are  somewhat  different  than  those  found  in  the  litera- 
ture primarily  due  to  their  dependence  on  both  radius  and  size  ratio. 

It  is  therefore  important  at  this  time  to  determine  the  effect  such  a 
difference  can  have  on  the  development  of  rain  in  a cloud.  For  this  purpose 
we  used  the  stochastic  numerical  model  of  Scott  and  Levin  (1975).  This 
model  treats  the  development  of  the  drop  spectra  with  57  size  classes  as  a 
function  of  time  in  an  infinite  cloud  model.  The  model  was  tested  with 
^coalesce  “ ^ compared  with  the  spectra  obtained  by  the  use  of  the 
coalescence  efficiency  as  described  by  Whelpdale  and  List  (1971) (Eq.  1-2) 
and  by  our  results  (Eq.  2-1).  The  results  of  the  size  spectra  after  600 
seconds  for  these  three  coalescence  efficiencies  are  shown  in  Fig.  2-14. 

The  initial  size  distribution  was  taken  as  a gamma  distribution  with  radius 
dispersion  of  0,28.  The  liquid  water  content  was  taken  as  3 gr/m^  and 
initial  concentration  of  100  drop/cm^.  In  the  graph  G represents  the  mass 
of  water  in  grams  per  m^  of  air  per  AJin  r . It  is  seen  that  for  E = 1 the 
spectrum  develops  very  rapidly.  When  Whelpdale  and  List's  (1971)  equation 
is  used  a considerable  decrease  in  the  conversion  of  cloud  water  to  rain 
water  is  observed.  When  the  present  results  are  used  an  intermediate  rate 
of  auto conversion  is  obtained.  These  results  are  very  interesting  since 
the  efficiency  obtained  here  is  lower  than  that  of  Whelpdale  and  List  at 
least  at  low  p ratios,  but  is  higher  at  other  p ratios.  The  fact  that  the 
present  results  lead  to  faster  growth  of  precipitation  elements  means  that 
interactions  of  drops  of  size  ratios  between  0,5-0. 9 are  the  most  important 
for  the  rain  development. 


CHAPTER  3:  THEORETICAL  CONSIDERATIONS 

3-1,  INTRODUCTION 

Unfortunately  there  is  no  theoretical  explanation  to  the  problem  of 
coalescence.  From  the  previous  work  in  the  literature  and  from  our  work 
described  in  Chapter  2 some  general  features  are  observed.  These  can  be 
classified  into  effects  or  forces  hindering  coalescence  and  those  that 
aid  coalescence.  In  the  following  few  pages  we  will  discuss  these  effects 
or  forces  in  order  to  justify  some  of  the  assumptions  which  will  be  made 
in  the  theoretical  treatment  that  follows. 

a.  Effects  that  hinder  coalescence 
i.  The  air  gap 

If  two  rigid  spheres  were  to  approach  each  other 
the  intervening  air  between  them  has  no  difficulty  in  escaping  and  clearing 
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the  way  for  actual  contact.  This  contact,  of  course,  will  occur  at  a 
single  point  on  the  sphere's  surface.  Water  drops,  however,  are  not  rigid 
and  as  forces  act  on  them  they  deform  to  conform  to  the  lowest  energy 
level.  When  two  water  drops  approach  each  other  the  pressure  in  the 
intervening  air  between  them  increases  and  forces  the  drops  to  deform. 

This  deformation  in  turn  hinders  further  drainage  of  the  air  and  slows 
down  the  approaching  drops.  About  the  existence  of  this  air  layer  one 
can  learn  from  our  work,  from  that  of  Whelpdale  and  List  (1971)  and  others 
in  which  air  bubbles  were  observed  inside  the  drop  after  coalescence. 

Also  the  works  of  Canals  (1960)  and  Selman  (1961)  reported  an  increase 
in  coalescence  with  a decrease  in  the  environmental  pressure.  The  force 
exerted  by  this  air  layer  is  the  viscous  force,  F , which  will  be  dealt 
with  in  greater  lengths  later. 

ii.  Surface  einergy 

Swinbank  (1947)  and  Browne  et  al  (1954)  suggested  that  if  a 
spherical  drop  will  deform  in  order  to  coalesce  with  another,  there  is 
a need  for  additional  surface  energy.  This  additional  surface  energy, 
even  though  small,  can  hinder  coalescence.  Similar  arguments  appear  in 
the  work  of  Cotton  and  Gokhale  (1967).  We  feel  that  the  deformation  does 
not  occur  for  the  coalescence  but  rather  because  of  it  as  will  be  seen 
later.  In  fact  the  surface  energy  acts  as  a driving  force  in  completing 
the  coalescence  once  contact  is  made. 

iii.  Aging  of  the  gas-liquid  interface 

The  time  lapse  between  the  formation  of  the  liquid  surface  and 
the  experiment  can  result  in  adsorption  of  gas  molecules  into  the  surface 
and  a decrease  in  the  ability  of  the  drops  to  coalesce  (eg,  Charles  and 
Mason,  1960,  Park,  1970).  Another  explanation  could  be  the  presence  of 
surfactants  on  the  drop's  surface  which  could  prevent  coalescence, 

iv.  Drop  oscillation  after  impact 

Experimental  observations  showed  that  after  impact  drops  behave 
as  if  they  Impacted  a solid  surface.  The  drops  often  get  a disk  like 
deformation  which  Increases  in  sizes,  reaches  a maximum  and  then  decreases 
As  the  drop  gets  rounded  it  acquires  a momentum  in  the  opposite  direction, 
just  like  an  elastic  ball  bouncing  off  a solid  wall.  Foote  (1975)  used  a 
numerical  method  to  calculate  the  drops'  shape  as  a function  of  time. 

From  his  computation  he  got  the  bouncing  time,  tg,  which  is  the  characteristic 
time  for  the  drop's  oscillation.  He  found  that  tg  is  very  similar  to  the 
time  of  the  first  mode  of  oscillation.  This  oscillation  is  refered  to  as 
the  oblate-prolate  mode  with  a period 
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where  a is  the  surface  tension,  R is  the  radius  of  the  drop  and  p its 
density.  It  therefore  means  that  after  from  the  initial  point  of 
collision  if  no  coalescence  occurs  the  drops  will  move  apart.  In  other 
words  according  to  Foote  the  impact  can  be  divided  into  two  processes. 

The  first  is  the  impact  with  the  solid  like  air  layer  and  the  second  is 
the  drainage  of  the  air  film.  This  drainage  competes  with  the  oscillation  of 
the  drops.  If  the  drainage  is  faster  than  coalescence  will  occur. 

Otherwise  a bounce  will  result. 

V.  The  electrical  double  layer 

Due  to  the  polar  nature  of  the  water  molecules  they  tend  to 
arrange  themselves  at  the  water  - air  interface  with  their  oxygen  atoms 
pointing  outward.  This  arrangement  creates  an  electrical  double  layer 
at  the  interface  which  acts  as  a negatively  charged  surface  at  small 
distances.  Hence  when  two  drops  approach  eaca  other  the  electrical  double 
layers  on  both  surfaces  repulse  each  other.  This  force  might  be  important 
at  very  small  distances  of  separation  and  will  be  discussed  later  on. 


b.  Effects  that  assist  coalescence 
1.  Electrical  effects 

Many  investigators  reported  on  the  effects  of  electrical 
charges  and  electric  fields  on  the  coalescence  of  drops.  Most  find 
that  below  a certain  charge,  bouncing  occurs  as  if  no  charge  was  present. 
However,  when  the  charge  exceeds  about  10~*  esu  the  coalescence  and  partial 
coalescence  probabilities  increase  (e.g.  Jayaratne  and  Mason,  1964,  Telford 
and  Thorndike,  1961,  Wheldpale  and  List,  1971  and  others). 

Electrical  forces  can  affect  coalescence  in  the  following  ways:  (1) 

Decrease  the  deformation  and  so  decrease  the  viscous  force,  Fy,  (2) 

Increase  (or  decrease)  the  attractive  forces.  (3)  Cause  elongation  of 

the  close  surfaces  and  there  by  increase  the  probability  for  initial  contact. 

li.  Oscillations  on  the  deformed  surfaces 

If  waves  or  disturbances  exist  on  the  approaching  deformed  surfaces 
they  can  assist  in  the  coalescence  process  by  decreasing  the  distance  of 
separation.  If  peaks  of  such  waves  are  opposite  each  other  they  can  actually 
start  the  coalescence  process  even  though  the  average  separation  distance  is 
large.  In  other  words  the  minimum  distance  of  separation  is  the  determining 
factor  in  the  coalescence  and  not  the  average  distance. 


iii.  Internal  ciruclation 

As  the  drops  fall  in  free  fall  in  the  atmosphere  an  internal 
circulation  within  them  is  developed  as  a result  of  the  stresses  exerted 
on  their  surfaces  by  the  air  flow.  This  internal  circulation  enhances 

the  drainage  of  the  air  film  between  the  drops  due  to  the  fact  that  the 
air  velocity  at  the  air-water  interface  within  this  region  does  . 
not  vanish.  Foote  (1975)  has  shown  that  coalescence  between  two  equal 
size  large  drops  occurred  only  when  internal  circulation  was  included. 

Iv.  Van  der  Waals  force 

This  attractive  force  could  act  between  the  approaching  drops 
at  very  small  separation  distances.  However,  calculations  show  that  for 
distances  larger  than  O.ly  it  could  be  neglected  (Whelpdale  and  List,  1971). 

c.  The  shape  of  the  deformation  and  its  size 

In  our  previous  discussion  we  mentioned  that  the  deformation  of  the 
approaching  surfaces  resembles  two  parallel  disks.  This  is,  of  course,  an 
approximation.  Allan  et  al  (1961),  Lindblad  (1964)  and  others  have  observed 
a dimple  like  deformation.  Burrill  and  Woods  (1969)  have  observed  two  such 
dimples.  Attempts  to  describe  this  deformation  as  a function  of  time  were 
also  made  (e.g.  Burrill  and  Woods,  1969).  Charles  and  Mason  (1960)  developed 
an  equation  for  the  viscous  force  under  any  symmetric  deformation. 

Despite  all  these  attempts  to  describe  the  forces  under  different 
deformations,  the  most  commonly  used  is  the  disk  like  deformation  (e.g. 
Plumlee,  1964,  Foote  1975).  The  justification  for  the  flat  deformation 
is  based  on  the  facts  that  some  experimental  observations  can  be  explained 
by  this  assumption  (MacKay  and  Mason,  1963) and  some  theoretical  analysis 
(Frankel  and  Mysels,  1961)  show  agreement  to  within  2%  with  experiments  when 
this  assumption  is  used.  In  the  present  work  the  disk-like  deformation  is 
assumed  especially  for  the  computations  of  the  viscous  force,  F^. 

In  addition  to  the  shape  of  the  deformation,  its  size  is  also  of  great 
Importance.  Derjaguin  and  Kussakov  (1939)  suggested  that  the  radius  of 
deformation  is  proportional  to  the  square  of  the  radius  of  the  drop  Itself 
and  is  a constant  throughout  the  Interaction.  However,  their  results  only 
apply  to  slowly  moving  drops. 

At  higher  velocities  the  deformation  is  a function  of  time  as  was  shown 
by  Plumlee  (1964)  and  by  Foote  (1975). 

In  the  present  work  we  used  two  models  to  describe  the  deformation.  One 
in  which  the  deformation  is  constant  during  the  interaction  and  the  other 
where  it  Increases  as  the  separation  distance  decreases. 
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d.  The  coalescence  distance 

The  coalescence  distance  is  defined  as  the  observed  minimum  distance 
between  the  drops'  surfaces  before  coalescence  occurs.  Observations  of  the 
coalescence  distance  are  very  difficult  to  obtain.  Llndbald  (1964)  observed 
it  to  be  less  than  1450  A while  Allan  et  al  (1961)  reported  that  the  coale- 
scence distance  was  about  900  A.  Whelpdale  and  List  (1971)  obtained  by  energy 
consideration  a coalescence  distance  of  between  0.1-ly.  Sartor  and  Abbott 
(1972)  computed  the  coalescence  distance  by  observation  of  the  air  bubble 
trapped  during  the  coalescence.  They  obtained  about  0.16y  for  the  coale- 
scence distance.  Based  on  all  these  observations  we  decided  to  use  a coale- 
scence distance  of  O.lji. 


3-2.  HEAD  ON  COALESCENCE 
a.  The  Model 

In  this  model  the  coalescence  of  two  liquid  drops  is  investigated. 

The  Model  is  primarily  aimed  at  describing  the  effect  of  the  different 
forces  on  the  approach  of  an  upward  moving  droplet  toward  a suspended 
drop,  as  was  discussed  in  the  experimental  section  of  the  report.  However, 
as  will  be  shown  later,  the  results  can  be  applicable  to  a variety  of 
other  type  of  experiments,  e.g.  interactions  at  relatively  high  impact 
velocities  such  as  Impact  of  two  drops  of  different  sizes  at  free  fall, 
impact  into  a liquid  surface  and  others.  It  is  not  applicable,  however, 
to  cases  where  a drop  is  placed  at  rest  on  a flat  liquid  surface  or  to 
cases  where  interactions  of  two  large  drops  occur  at  free  fall. 

The  model  assumes  that  the  origin  of  the  coordinate  system  rests  at 
the  center  of  the  large  drop  and  that  it  rests  with  respect  to  the  surrounding 
air.  The  smaller  drop  moves  vertically  upward  at  a velocity  AV  which  is 
the  difterence  in  th^"%«:minal  velocities  of  the  drops.  Since  the  model 
describes  the  approach  of  the  two  drops  it  is  possible  to  determine  whether 
coalescence  occurs.  In  the  next  chapter  the  same  method  will  be  used  to 
determine  conditions  for  coalescence  in  non-head-on  interactions.  From 
all  these  coalescence  efficiencies  for  different  drop  sizes  will  be  deter- 
mined. In  order  to  develop  the  equation  of  motion  we  first  have  to  find  an 
expression  describing  the  forces  acting  on  the  drops.  The  viscous  force, 
the  force  of  the  electrical  double  layer.  Van  der  Waais  force  ana  electrical 
forces  all  depend  on  the  shape  of  the  drops.  Hence,  the  approach  of  the 
drops  also  depends  on  their  shape. 

To  our  knowledge  there  is  no  way  in  which  the  shape  of  approaching 
liquid  drops  can  be  expressed  mathematically.  All  the  expressions  found 
deal  with  a quasi  stationary  problem  in  which  the  drops  are  almost  at  rest 
with  respect  to  each  other.  For  this  reason  we  chose  two  simple  deformations 
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to  express  the  changes  that  take  place  on  the  surfaces  the  approaching 
drops.  However,  the  equations  we  will  propose  can  be  easily  changed  when 
a better  expression  for  the  deformation  is  found. 

b.  The  Forces  and  the  deformation.  The  forces  acting  on  a droplet 

which  moves  upward  toward  the  center  of  a larger  suspended  drop  are 

i.  The  viscous  force  F^.  This  force  arises  as  a result  of  the 
resistance  of  air  to  the  motion  of  the  drop.  At  large  distances  of 
separation,  this  force  depends  on  the  size  of  the  smaller  drop,  on  its 
velocity,  and  on  its  instantaneous  deformation.  At  small  distances  of 
separation,  F^  depends  on  the  deformation  of  the  two  drops  and  on  the 
distance  between  them. 

ii.  The  gravitation  force,  F . This  force  acts  in  a direction  so 
as  to  prevent  the  upward  motion  of  fhe  droplet. 

iii.  The  electrical  double  layer  force,  F This  repulsive  force 
results  from  the  electrical  dipoles  which  align  themselves  along  the 
water-air  interface.  They  are  arranged  so  that  at  small  distances  of 
separation  the  water  surfaces  act  as  two  negatively  charged  surfaces. 

This  force  depends  on  the  size  of  the  two  drops  and  their  distance  of 
separation. 


iv.  The  London-Van-der-Waals  force,  This  attractive  force 

depends  on  the  deformation  of  the  drops  and  the  distance  between  them. 

Let  D be  the  Instantaneous  minimum  distance  between  the  two 
drops.  We  are  interested  in  effects  which  take  place  at  distances  larger 
than  10"®  cm  since  it  is  widely  reported  that  coalescence  occurs  at 
distances  of  the  order  of  lO'*  cm  (the  mean  free  path  of  air  molecules 
at  room  temperatures).  According  to  Adamson  (1967),  F^^  and  Fj^y  per  unit 
area  are  given  by 

Fdi  = (8.9*10  ’/D^)  (dynes/cm^)  (3-2) 

Fj^y  - (10“^  ^/6irD^)  (dynes/cm^)  (3-3) 

so  that 

- 54tt10®D  (3-4) 

which  Implies  that  F^^^  and  Fj^y  are  equal  when  D = 10"®/54'TT  = 5.9' 10~®  cm 
so  that  for  any  D larger  than  10"®  cm,  Fj^y  <<  Fdl*  that  Fj^y  can  be 
neglected. 
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The  three  forces  F„,  F and  F^^  are  all  repulsive  forces  in  the  present 
situation.  Qualitatively,  the  behavior  of  the  small  drop  and  the  forces  on 
it  can  be  described  as  follows:  the  small  drop  moves  upward  it  an  initial 
velocity  Vq  which  decreases  as  it  approaches  the  larger  drop.  Where  no 
coalescence  occurs  the  small  drop  stops,  the  pressure  between  the  drops 
decreases  and  the  deformation  of  the  drops  decreases  as  well.  This  decrease 
in  the  viscous  force  may  result  in  oscillations  about  the  equilibrium 
spherical  shape  of  the  drops  which  may  enhance  the  probability  for  coalescence. 
On  the  other  hand,  since  at  that  point  Fg  and  F^j^  still  act,  the  drops  are 
forced  apart.  These  last  two  effects  compete  with  each  other  and,  hence, 
there  is  a need  to  treat  them  both  simultaneously.  In  this  report  we  deal 
only  with  the  motion  of  the  small  drop  to  the  point  at  which  it  stops  and 
neglect  the  effect  of  the  oscillations. 


The  equation  of  motion  of  the  small  drop  can  be  written  as: 

= Fdl  "g  "V 

Since  F^  and  F strongly  depend  on  the  shape  of  the  approaching  surfaces, 
we  first  have  to  discuss  the  deformation  of  the  drops. 

V.  The  Deformation.  The  deformation  of  the  drop  depends  on  its 
size,  on  the  distribution  of  pressures  around  it,  and  on  its  surface 
tension.  It  is  an  impossible  task  to  evaluate  the  pressure  distribution 
of  the,  approaching  drops  at  each  point.  As  a reasonable  approximation 
it  is  usually  assumed  that  the  approaching  surfaces  are  planes  (even  though 
it  is  believed  that  the  larger  drops  tend  to  get  a dimple-like  deformation), 
Frankel  and  I^sels  (1962)  showed  that  the  approximation  of  a planar  deforma- 
tion of  a drop  approaching  a water  surface  gives  reasonable  results  to 
within  2%  error. 

We  assume  a planar  deformation  proportional  to  the  radius  of  the 
smaller  drop.  Two  deformation  are  attempted: 

I.  Rp  = aRg  (3-5) 

II.  Rp  « 6Rg/D  (3-6) 

where  D is  the  distance  between  the  deformed  surfaces,  R^  is  the  radius  of 
deformation  and  R^  is  the  radius  of  the  small  drop,  a and  B are  constants. 
Since  the  deformation  of  the  large  drop  is  a direct  result  of  the  approach 
of  the  small  drop,  we  assume  that  R^  is  the  same  for  both  drops. 

vi.  The  Viscous  Force.  This  force  acts  mainly  due  to  the  fact  that 
the  deformed  surfaces  approach  each  other.  We  can  therefore  use  the  force 
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between  two  planes  of  radius  Rq  moving  toward  each  other  at  a velocity  V 
(see  Landau  and  Lifshltz,  1959): 


BTtnRp  ^ 

2 D3 


(3-7) 


It  is  interesting  to  note  that  this  was  the  starting  point  for  many 
previous  ^'orks  (e.g.  Charles  and  Mason,  1960).  However,  they  assumed 
that  Fy  in  addition  to  Rp  are  constant  and  so  arrived  at  the  Stef an-Reynolds 
equation : 


- "o  = 


3TrnRp“ 
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However  the  assumption  of  a constant  viscous  force  is  unreasonable  as 
will  be  shown  later.  In  the  present  work,  contrary  to  most  other  works, 
we  are  not  assuming  the  quaSl-stationary  conditions  and  the  equations  are 
solved  without  neglecting  the  inertial  force  mD.  Foote  (1971)  has  also 
considered  inertial  forces  but  only  with  respect  to  the  plane  of  deformation 
and  not  with  respect  to  the  other  drop. 


vii.  The  Electrical  Double  Layer  Forces.  When  a planar  deformation 
is  considered  it  is  reasonable  to  assume  that  only  the  deformed  portion  of 
the  drop  contributes  to  the  electrical  double  layer  force,  and  according 
to  Adamson  (1967)  it  can  be  expressed  as  follows: 


8.9  • 10-7t:R^2 


dynes 


(3-8) 


viii.  The  force  of  gravity.  If  the  radius  of  the  small  droplet  is 
the  force  of  gravity  on  it  would  be 

Fg»-|TTRjg  (3-9) 

This  force  can  enhance  coalescence  or  inhibit  it  depending  on  the  interaction. 
In  the  case  where  a small  droplet  impacts  the  lower  hemisphere  of  a large 
drop,  the  force  of  gravity  tends  to  oppose  coalescence.  In  impact  of  drops 
onto  flat  liquid  surfaces  it  assists  in  the  coalescence.  In  all  cases  tested 
here  where  the  impact  speed  was  comparable  to  the  relative  terminal  velocities 
of  drops  between  20-2000vi  the  effect  of  gravity  was  very  small  and  could  be 
neglected . 


c.  The  equation  of  motion 

The  model  deals  with  a small  drop  moving  upward  toward  a large 
drop.  The  forces  Fy,  F^^^  and  F all  act  downward.  It  can  be  easily 
shown  that  the  connection  betwefn  r and  D can  be  written  as: 

r = D + /r2  - + /Rg  - (3-10) 


where  r is  the  distance  between  the  centers  of  the  drops  and  R^  is  the 
radius  of  the  large  drop.  The  equation  of  motion  of  the  small  droplet 
can  now  be  written  as 

m ^ = ZF  (3-11) 

s dt^ 

whe  re 


ZF 


F + F + F + F 
V DL  LV  g 


2D^  1)2 


10~^^ 

6D^ 


+ 


4TTr4 

3 


g 


(3-12) 


and  where  n is  the  viscosity  of  air,  D is  the  distance  between  the  approach 
ing  surfaces  and  m^  is  the  mass  of  the  small  droplet. 


It  is  worthwhile  at  this  stage  to  change  the  dependent  variable  from  r 
to  D.  By  the  use  of  Eq.  (3-10)  we  can  show 


d^ r _ 2 d^ r ^ d^  D dr 

dt^  dD^  dt^  dD 


(3-13) 


where 


v 


= 

dt 


is  the  relative  velocity  of  the  deformed  surfaces. 
If  we  use  the  first  deformation  criterion: 


R„  = aR 
D s 

we  get 

ill  = ili 

dt^  dt^ 


so  that  the  equation  of  motion  becomes 
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(3-14) 


If  the  second  defoimation  is  used  namely 

„ SRs 


'D  - D 


than  we  get 


- 1 + ® ( 

dD  D3  '■ 


(3-15) 


dll » _ ( ..  1.  + _JL_ ) _ !A( 1 H 


__i + 1 ) (3-16) 


and  the  equation  of  motion  looks  like 


ilD,_Lj.  (dD  D).y2d?r  /dr  (3_17) 

dt2  mg  ^ dt/  ^ dD2  V dD 


substitution  of  Eqs.  (3-15)  and  (3-16)  into  (3-17)  gives  the  equation  of 
motion  without  r. 


The  order  of  the  above  differential  equations  can  be  reduced  by  the 
use  of 


d D dv 
dt^  ° dD 


Using  this  relationship  we  get  for  R^ 


BRg 

and  for 


4^  = — E F(v,D) 
dD  m V * 

s 


(3-18) 


(^zr  (v.W  - v|^)/f 

O ' 


(3-19) 


This  equation  is  applicable  for  v ^ o 


Both  equations  were  solved  numerically  by  the  use  of  a forth  order 
Runge-Kutta  method.  For  small  values  of  v it  was  required  to  substantially 
reduce  the  step  size  in  order  to  achieve  better  accuracy.  A great 


Jlgt 


iA 
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improvement  was  obtained  by  subsituting 


( — 

^ dt  ’ 


(3-20) 


With  this  substitution  we  obtain: 


for  = aR 
D s 


^ = ZF(u.D) 


(3-21) 


R - S 


L 


^ - 0:f(u.d)  - 2 


d r j dr 


d^  r.,  /dr 


(3-22) 


where  both  — — and  — can  be  substituted  from  Eqs.  (3-15)  and  (3-16), 
dD^  do 

u has  the  meaning  of  kinetic  energy  when  R^^  = aR  is  used,  but  it 
has  no  simple  meaning  in  the  second  deformation  criterion. 

A solution  of  Eqs.  (3-18)  and  (3-19)  results  in  a negative  v 
which  implies  an  approach  of  the  drops.  At  a very  small  distance  of  separa- 
tion, D,  V may  get  positive  values,  but  the  viscous  force  is  not  defined 
for  such  a situation  hence  this  solution  is  not  physical,  u on  the  other 
hand  describes  the  relative  velocity  of  the  two  approaching  drops,  hence 
for  all  real  physical  situations  it  is  positive. 

d.  The  numerical  solution 

The  numerical  solution  to  the  approach  velocity,  v,  as  a function 
of  separation  distance  for  different  drops  with  the  same  initial  velocity 
(-600  cm/sec)  is  shown  in  Fig.  3-1.  Part  (a)  describes  the  velocity  for 
deformation  Rq  = 0.5  Rg  while  part  (b)  is  for  Rp  = 0.5  x 10“^Rg/D,  The 
following  observations  can  be  made: 

(1)  All  the  curves  decrease  monotonically 

(2)  At  large  separation  distance  the  velocities  are  all  the 
same  and  are  equal  to  the  initial  velocity,  V^, 

(3)  Changes  in  the  initial  separation  distance,  do  not 
alter  the  results  as  long  as  is  large. 

(4)  The  larger  is  the  small  drop  the  more  rapidly  will  its 
approach  velocity  decrease.. 
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(5)  Each  curve  has  a certain  Dg  to  which  it  approaches.  This 
distance  is  the  stop  distance  of  the  deformed  surfaces, 

(6)  The  larger  is  the  smaller  drop  the  larger  is  Dg, 

(7)  The  slopesof  the  curves  increase  as  is  approached. 

Fig.  3-lb  represents  the  same  characteristics  as  above 
except  that  the  slope  increases  to  a distance  D > Dg  and 
decreases  as  Dg  is  further  approached. 

As  [V  I increases  so  does  v at  every  D.  This  results  in  a decrease 
in  the  value  of  Dg.  This  is  demonstrated  in  Fig.  3-2  for  the  case  - 50ym. 
Again  Fig.  3-2a  represents  the  case  Rq  ■ fig*  3-2b  represents 

the  deformation  Rq  ••  0.5  x 10"^Rg/D. 

The  head-on  approach  depends  on  Rg  and  and  not  on  Rj^.  However, 
since  we  are  interested  in  interactions  of  drops  falling  at  terminal 
velocities,  R^  and  Rg  determine  and  therefore  influence  the  approach. 

a and  B,  the  coefficients  of  the  deformation  radius,  strongly  Influence 
the  approach  of  the  drops.  As  a and  B get  smaller  so  do  the  deformations 
and  the  velocity  at  each  D is  larger  and  results  in  closer  approach  of 
the  drops  and  a smaller  Dg.  This  is  demonstrated  in  Fig.  3-3  for  the  case 
Rg  ■ 50y,  R,  ••  250y  and  their  relative  terminal  velocities.  Part  (a) 
represents  the  first  deformation  criterion  R ■ aRg  where  a varies  between 
0.3  to  0.7.  As  can  be  seen  a has  a very  strong  effect  on  the  approach  of 
the  drops.  Part  (b)  represents  the  same  situation  for  a deformation 
Rp  - 6Rg/D  where  B varies  between  3 x 10“®  to  7 x 10”®.  The  effect  of  B 
on  the  approach  is  smaller  and  less  drastic  than  that  of  a. 

e.  A comparison  between  the  forces 

From  the  numerical  computations,  values  of  F^  can  be  obtained 
and  compared  to  the  other  forces.  Fig.  3-4  demonstrate  this  comparison 
for  the  case  where  R - 50y,  R^  " 250y,V^  is  their  relative  terminal 
velocities  and  Rp  =•  0.5  Rg.  It  is  observed  that  the  maximum  value  of  Fy 
is  larger  by  a few  orders  of  magnitude  from  the  values  of  all  the  other 
forces  throughout  the  important  range. 

We  define  Dj.j  to  be  the  point  at  which  F„  is  a maximum.  For  D > 10  Dj^ 
the  velocity  does  not  change  and  we  are  justified  in  assuming  mr  “ o. 

In  the  range  Dj^  < D < 10D[.(  F^  is  larger  than  all  the  other  forces  and  so 
they  can  be  neglected. 

For  D around  D^^  the  viscous  force  is  much  larger  than  all  the  other 
forces.  As  D gets  very  close  to  D^  the  other  forces  can  become  Important 
and  there  is  a need  to  consider  them  as  well.  However,  this  range  is  very 


D,  ClISTArJCE  OF  SEPARATION  CCIVlD 

Fig.  3-2a.  The  dependence  of  approach  velocity  on  the  distance  of  separation  for 
different  initial  velocities.  Here  R-  = 3.5R  , R * 50y  and  R,  » ZSOu. 
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close  to  Dg  and  therefore  very  small, 
to  a first  approximation  for  all  values 
has  to  be  considered. 


It,  therefore,  can  be  concluded  that 
of  D > Dg  only  the  viscous  force 


Is  <tefor«atlo„  criterion 

»sJs-»canl:!  " the^custlon  of  t„,tlo„  L be 


f.  The  analytic  solution 

The  eq 
can  be  written; 


ThCeequation  of  motion  with  the  assumption  of  deformation 


‘ OiR„ 


m d D 
dt^ 


1 d_D 
2 r>3  dt 


(3-23) 


and  can  be  rewritten  as 


where 


37m  (aRg)** 


(3-24) 

(3-25) 


and  where  a dot  represents  a time  derivative.  By  writing 


D « V 


dv 

dD 


we  get 

dv  _ A 
dD  ■ jj3 

3y  Integration  we  get 


If  we  define 


W - V - 
1 o 


-A-  + 

21^  2D^ 

(3-26) 

A 

2^0 

(3-2  7) 

t 

► 

! ‘ » 


I i 

I 


- 55  - 


we  can  re-write  Eq.  (3-26)  as 

A 


W.  + 


(3-2  8) 


Since  we  deal  with  approaching  surfaces  the  equation  had  s meaning  only 
if 


D > /■ 

2Wj 

has  the  meaning  of  the  velocity  that  the  small  drop  was  to  have  at 
infinity  in  order  to  reach  point  at  a velocity  V^i, 

The  analytic  solution  provides  a better  understanding  of  the  approach 
problem,  for  example: 


(1)  If  D » 
o 


2v„ 


then 


(3-29) 


^ = '^o  2D^ 


Here  D does  not  appear  at  all  and  therefore  it  implies  that  in  this  case 
the  approach  does  not  depend  on  D^, 


/ A 

(2)  For  D >>  /-  -r — the  velocity  is  almost  const£int  v « v 

4 V C 

O 

(3)  The  value  of  can  be  found.  We  know  that 

lim  V = o 
D -»■  D 


Therefore 


2Wj 


(3-30) 


If  we  now  define  ■*  D/Dg  and  = v/w^  we  get  a relation  between  a 
normalized  velocity  and  normalized  distance 


V - 1 - ^ 

Fig.  3-5  represents  Eq.  (3-31), 


(3-31) 


i_.  _ Jlgifc, 


Fig.  3-S.  KonaalizeU  velocity  and  viscous  force  as  a function  of  the  nornalized 
distance  of  separation. 
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g.  The  viscous  force 

Substitution  of  Eqs.  (3-25)  and  (3-28)  into  Eq.  (3-6)  gives 


F,(D)  . (-  ^ ^ 


It  is  therefore  easy  to  see  that 


W ■ 


m AW.  . , 

^ i / I J_  s 

^ iTS’  D,^  ^ 
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N 


(3-32) 


m i^Wf 

( ""  '*•  ^ 

is  the  viscous  force  that  would  have 


The  coefficient 

acted  between  the  drops  at  a distance  if  the  drop  did  not  lose  part 
of  its  speed.  If  we  define  a normalized  viscous  force 


(- 


dJ 


(3-33) 


we  get 


F = -i i- 

N N 


(3-34) 


This  equation  is  also  represented  in  Fig,  3-5.  F becomes  maximum  at 


/|  D„  = 1.; 


3 “s  -291 


At  that  point  F^  is: 


=-2(3  )3/2 
max  5 '■  5 ^ ^3 


when  F is  a maximum  V = 2/5.  As  the  two  drops  approach  each  other  and 
get  ve^  close  to  D , the  viscous  force  decreases  rapidly. 


h.  The  approach  time 

It  is  sometimes  interesting  to  calculate  the  time  it  takes  fo/ 
the  droplet  to  approach  to  within  a distance  D from  another  drop  under  the 
action  of  its  intertial  and  viscous  forces. 
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The  time  can  be  found  by: 


Substitution  of  V from  Eq.  (3-28)  and  Integration  results: 


D-d  ° Ds 


(3-35) 


where  T Is  a time  constant  defined  as: 


2(-W^) 


According  to  this  equation  the  drop  gets  closer  to  D but  reaches  It  only 
after  Infinite  time.  However,  one  has  to  remember  tEat  the  analytic  solu- 
tion Is  not  applicable  at  distances  close  to  Dg,  since  at  this  distance 
other  forces  have  to  be  considered. 

1.  The  range  In  which  the  analytic  solution  Is  valid 

(a)  For  the  analytic  solution  to  be  physically  valid  Fy  » F 
and  also  Fy  » Fdl  (D^) • For  these  conditions  we  find 


i 5.4  /r 
1 s 


- W.  i 0.0103 

^ % 


(3-36) 


respectively.  For  the  range  20y  ^ Rg  < 2000y  It  Is  enough  that  > 5.15 
cm/sec  for  these  conditions  to  be  fulfilled. 

(b)  Lets  assume  that  there  Is  a range  D^]  In  which  the  analytic 
solution  Is  valid  and  lets  try  and  find  It. 

For  large  D the  velocity  v “ Vq  If  within  the  range  D to  D^^  there  was 
no  appreciable  lose  of  velocity.  For  large  D Fdl  « Fg  and  hence  D^  will 
be  determined  by: 

D - 1.291  D 

g ( ° , ^ 100 
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From  all  these  we  can  conclude  that  the  analytic  solution  Is  valid  for  most 
cases  considered  where  the  Impact  velocity  Is  higher  than  51  cm/sec  and  Is 
fairly  good  for  even  lower  values  of  However,  for  lower  Impact  speeds 

and  Impact  at  angles  different  from  head-on  a numerical  solution  Is  needed. 


CHAPTER  4:  NON -HEAD-ON  APPROACH  OF  TWO  DROPS 


In  this  chapter  we  follow  the  motion  of  two  drops  as  they  approach  In 
a non-head  on  collision  course.  From  considerations  of  this  approach  It 
Is  possible  to  draw  conclusions  regarding  the  coalescence  Itself. 

a.  The  equations  of  motion 

The  origin  of  the  coordinate  system  was  centered  at  the  center  of 
mass  of  the  large  drop  with  the  assumption  that  the  coordinate  system  Is 
Inertial.  The  equations  of  motion  of  the  small  drop  In  polar  coordinates 
are : 


(4-1) 

m r - m r0^  - F„ 
s s R 

(4-2) 

where  r Is  the  distance  between  the  centers  of  the  two  masses,  m the  mass 
of  the  smaller  drop  and  F^  and  Fj^  are  the  tangential  and  radial  forces, 
respectively.  In  our  calculations  here  we  assume;  (1)  That  the  planes  of 
deformations  on  the  two  drops  are  parallel.  This  assumption  causes  the 
forces  Fy,  Fp^  and  Fpy  to  act  radially.  In  order  to  get  Fy  we  substitute 
In  Eq.  3-12  the  relative  velocity  of  the  deformation  planes,  v = dD/dt  in 
the  radial  direction.  (2)  Since  in  the  tangential  direction  the  drop  looks 
spherical  » F^  and  hence  Fj.  will  be  neglected. 

Here  we  have  to  define  a few  terms: 

Impact  velocity  - is  the  relative  velocity  of  the  drops  before  impact. 

Impact  trajectory  - is  the  line  drawn  before  impact  from  the  center  of 
the  small  drop  In  the  direction  of  motion. 

Impact  point  A - The  Intersection  of  the  impact  trajectory  with  a sphere 
of  radius  equal  to  the  sum  of  the  radii  of  the  drops  and  centered  on  the  center 
of  the  large  drop  (see  Fig.  4-1) . 

Impact  angle  9.  - is  the  angle  between  the  direction  of  the  impact  velocity 
of  the  small  drop  and  the  line  connecting  the  impact  point  with  the  center  of 
the  large  drop. 


I. 


Impact  point  according  to  Whelpdale  and  List  (1971)  - is  the  inter- 
section of  the  impact  trajectory  of  the  small  drop  with  the  large  drop. 
This  point  is  marked  as  B on  Fig.  4-1.  The  impact  angle  resulting  from 
this  last  definition  is  called  a^. 

Impact  distance  - is  the  distance  between  the  impact  trajectories 
of  the  drops.  The  impact  angle  and  the  impact  distance  provide  means  to 
express  how  close  to  head-on  is  the  impact 


^ ^s^  ®i 

x^  » R sin 


(4-3) 


Instantaneous  impact  distance,  x - is  the  distance  between  the  lines 
drawn  from  the  drops'  centers  in  the  direction  of  motion. 

As  long  as  the  drops  are  far  apart  x = x^  and  V » V^.  Effects  of 
collision  or  coalescence  cause  changes  in  x and  V. 

From  Ea . (4-1)  we  get 

r^  0 = C 

where  C is  a constant  of  the  motion.  For  V^(the  tangential  component  of 
V)  we  have: 


Vj.  = r 0 

so  that  at  any  time 
C = r 

Vt 

and  tan  a = — and  tana=  x/r. 

Therefore  at  any  instant: 

C = Vx 

If  the  initial  conditions  are  marked  with  subscript  zero: 

C=  V X (4-4) 

o o 

The  initial  conditions  have  to  be  chosen  at  such  distances  that  the  effects 
of  collisions  are  small.  For  impact  on  a stationary  target  it  is  often 
assumed  thatitlstrue  at  any  distance,  and  the  initial  point  can  be  taken  as 
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X = X.  and  V = V . . Under  these  conditions 
o i o i 


(4-5) 


we  now  get 


m 


r 02  = 


m 


m (V  X 
S 0 0 


S 


and  substituting  into  Eq.  4-2  we  get 


m 


m (V  X )' 

S 0 0 


(4-6) 


where 


p>  j.  F + F 

'v  ^ 'dl  'lv 


The  gravitational  force,  F , is  not  necessarily  in  the  radial  direction, 
however  it  can  be  generally  neglected. 

Eqs . (4-])  and  (4-2)  can  be  divided  to  a tangential  equation  whose  solu- 
tion is : 

Vt(r)  - ^ (4-7) 


and  to  a radial  equation  (Eq.  4-6)  which  can  be  solved  numerically. 

The  second  term  on  the  right  hand  side  of  Eq.  (4-6)  is  actually  the 
centrifugal  force,  F^,  which  inhibits  the  approach  of  the  drops.  F^, 
increases  as  Xq  gets  larger.  For  Xq  = o the  equation  becomes  identical 
to  a head-on  collision. 


b.  The  numerical  solution  of  the  radial  equation 


Figs.  (4-2)  and  (4-3)  present  the  variation  of  v,  the  relative 
velocity  of  the  deformation  surfaces,  as  a function  of  D,  where  the 
Impact  distance  is  taken  as  a parameter.  The  graphs  were  drawn  for 
Rj)  ■ 0.5  Rg  and  Rp  * 0.5  x 10“®  Rg/D  respectively,  and  for  Rg  = 50y, 

Rj^  * 250y  and  as  the  relative  terminal  velocities  of  these  drops. 
Initially  (high  values  of  D)  the  approach  is  not  affected  by  the  presence 
of  Fy,  Fdl  and  F^y  So  that  the  approach  is  describes  by: 


m r 
8 


2 


Fig.  A— 3.  The  same  as  Fig.  4-2  except 
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The  solution  to  the  above  equation  Is: 


— = V ( 1 - ( 
dt  i'*  ^ r 


(4-8) 


When 


Rjj  • ^ “ dr/dt  one  gets 


v(D)  = ( 1 - 


)* 

[D  + (rJ  - -)-  (rJ  - 


(4-9) 


If  the  approach  were  to  continue  according  to  this  equation  then  for 
small  D (D  <<  (r5  - R^)^),  v would  not  depend  on  D but  would  approach  the 
horizontal  line  In  Fig.  4-2. 

BRfl 

For  Rjj  = — , Eq.  (4-8)  becomes 


v(D)  = 


(1  - 


(dr/dD) 


[D+(R^-r2)^+  (R^-R^)*^]^ 


(4-10) 


In  Eq.  (4-10)  a term  dr/dD  appears  and  is  given  in  Eq.  (3-15). 

For  large  values  of  D,  dr/dD  % 1 and  the  effect  of  the  radius  of 
deformation,  Rp,  is  small.  However,  as  the  drops  approach  each  other, 
the  other  forces,  and  especially  Fy,  dominate  and  reduce  the  value  of 
lv|.  The  larger  is  the  impact  distance  the  smaller  is  |v|  (for  the 
same  D) , and  the  larger  Dg(the  stop  distance)  will  be  when  the  drops  stop 
their  approach  (v  = o) . This  fact  is  true  for  both  deformations  even 
though  it  might  be  difficult  to  see  it  in  Fig.  4-3.  At  the  initial  point 
Eq.  (4-8)  looks  like: 


t 

i. 


t * 


when  dr/dt)  is  the  radial  component  of  the  initial  relative  velocity  of 
the  centers  of  the  drops. 

It  is  easy  to  see  that  if  > 7.1  then  '^)q  > 0.99  and  dr/dt) 
does  not  depend  on  x^^  or  on  the  exact  value  of  r . This  is  sufficient 
conditions  for  the  motion  to  be  independent  of  the  starting  place,  however, 
it  is  not  a necessary  condition.  For  smaller  values  of  r^, dr/dt)  will 
depend  on  r^,  but  as  long  as  Eq.  (4-8)  is  fulfilled  the  approach  will  still 
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be  independent  of  r.  In  fact  numerical  tests  indicate  that  the  approach 
of  the  drops  is  not  affected  by  the  initial  separation,  r^,  as  long  as 
Fv(ro)  is  small. 

Fig.  4-4  represents  the  radial  force  Fj^(»Fy)  and  the  centrifugal 
force,  Fj,,  as  a function  of  distance,  D,  for  Rg  ■ 50u  and  • 250y 
interacting  at  their  relative  terminal  velocity  V^. 


The  deformation  was  taken  as  Rd  - 0.5 
45' 


Rg  and  the  impact  angle 


The  other  force  Fg,  Fj^y  and  Fql  are  much  smaller  and  could 


% 

not  be  drawn  on  this  scale.  The  curves  indicate  that  down  to  a distance, 

D,  of  a few  microns  F dominates.  However,  at  smaller  separation  distances 
the  viscous  force  dominates.  Since  Fg  is  small  and  can  be  neglected 
(under  the  present  case  of  relatively  high  velocities)  the  solutions  can 
be  applied  to  different  interactions  irrespective  of  their  impact  direction. 


The  approach  of  the  drops  in  head  on  collision  was  independent  of 
the  size  of  the  larger  drop,  Rj^,  except  when  was  computed.  In  the 
general  case  R^  affects  the  centrifugal  force  through  V^: 


[D  + 


(R?  - 


(Ro  - 


Rn)r 


(4-11) 


It  indicates  that  as  Rj^  decreases  the  radial  velocity  will  vanish 
when  the  drops  are  further  apart. 


c.  Summary 

The  relative  velocity  of  the  deformation  planes  of  the  drops 
can  be  divided  into  tangential  and  radial  components.  The  tangential 
velocity  does  not  change  appreciably  during  the  approach  while  the 
radial  velocity  decreases  sharply.  If  no  coalescence  occurs  when 
V “ o,  the  small  drop  will  move  sideway  at  a velocity  equal  to  its 
tangential  velocity  at  that  point. 


I 

i 


The  larger  is  the  smaller  is  v and  the  lower  is  the  probability 
that  the  interaction  will  result  in  coalescence. 


DISTANCE  OF  SEPARATION,  O (CM) 


Fig.  4-A.  The  viscous  and  the  centrifugal  forces  as  a function  of  the  distance 
of  separation.  Here  Rg  - 50y,  Rl  ■ 250y,  Rp  - 0.5Rg  and  the  initial 
velocity  was  taken  as  the  relative  terminal  velocity. 
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CHAPTER  5.  THE  COALESCENCE  PROBLEM  AND  THE  COALESCENCE  EFFICIENCY 


The  problem  of  coalescence  concerns  with  the  question  of  whether 
two  colliding  drops  will  coalesce,  even  temporarily,  or  will  bounce 
off  each  other.  This  chapter  offers  a criterion  for  answering  this 
question,  for  determining  the  regions  of  coalescence  and  bouncing  from 
the  approach  of  two  drops,  for  determining  critical  contact  angles  and 
for  determining  impact  distances  and  coalescence  efficiencies. 


5-1.  An  answer  to  the  coalescence  problem 

Experimental  and  theoretical  considerations  suggest  that  the 
maximum  separation  distance  needed  before  coalescence  can  take  place 
is  about  O.ly.  We  therefore  take  this  distance  as  the  necessary  distance 
to  which  the  drops  have  to  approach.  If  the  approach  velocity  of  the 
two  deformed  surfaces,  v,  vanishes  at  a distance  larger  than  O.ly, 
bouncing  results.  If,  on  the  other  hand,  v has  a value  greater  than 
zero  at  = O.ly  the  interaction  results  in  coalescence. 

According  to  the  above  criterion  a computer  program  was  developed 
into  which  the  initial  conditions  of  the  motion  and  the  drops'  sizes 
were  fed  and  which  provided  the  answer  to  the  coalescence  problem  for 
each  pair  of  drops. 


5-2.  Regions  of  bouncing  and  coalescence 

Since  the  radial  component  of  the  velocity  is  smaller  for  larger 
impact  angles,  the  distance,  D^,  at  which  v = o is  larger.  This  implies 
that  head-on  collisions  that  result  in  separation  will  not  coalesce  at 
larger  inpact  angles.  These  cases  can  be  studied  from  the  results  of 
Chapter  3. 


The  more  complicated  problem  is  the  determination  of  the  angle  below 
which  coalescence  occurs  and  above  which  bouncing  takes  place. 

i.  Regions  of  bouncing  and  coalescence  in  (Rg,  |v^()  plane 
In  the  analytic  case  we  obtained 


T a"( 


-w. 


when  is  given  in  Eq.  (3-27).  For  known  a and  ri  and  large  D we  can 
write 

3 .2,  ^'^S  .4 


°s  = Ds(Rs*Vi)  = 4 “ 


Fig.  5-1  describes  schematically  Dg  versus  R and  |Vj^l  and  the  inter- 
cept of  the  function  with  the  plane  Dg  - O.lp.  This  intercept  is  a straight 
line  which  can  be  described  in  the  - Rg  plane  as 


|V^|  = (1.78  X 10^°nci‘*)Rg 


(5-1) 


Fig.  5-?  describes  this  equation  for  a few  values  of  Ot.  For  each  a there 
is  a curve  which  divides  the  plane  into  two  regions.  Above  the  curve  all 
interactions  result  in  coalescence.  Below  it  all  interactions  result  in 
separation . 

For  the  cases  where  an  analytic  solution  cannot  be  obtained,  a computer 
program  was  developed  which  provides  the  appropriate  curves  for  different  a. 

For  the  deformation  R^  = aRg  a full  correspondence  with  Fig,  5-2  was 
obtained.  For  Rp  = curves  presented  in  Fig.  5—3  where  obtained 

when  Rp  = lOO'Oy,  With  both  deformations  straight  lines  near  the  origin 
where  obtained  and  their  slope  was  found  to  increase  monotonically  with 
the  deformation  constant. 

For  every  | | the  point  (Iv^IjR^)  on  the  curve  gives  the  critical 
Rg  above  which  coalescence  takes  place  and  below  which  bouncing  occurs. 

In  the  same  way  for  each  Rg  there  is  a |V^1  above  which  coalescence  takes 
place.  This  result  explains  the  experimental  fact  that  in  head-on  collisions 
at  low  velocities,  coalescence  occurs  (Rayleigh,  1882,  Park,  1970  and 
Nelson  and  Gokhale,  1973). 

ii.  Regions  of  bouncing  and  coalescence  in  the  (pjRp)  plane 

If  one  is  interested  in  initial  velocities  which  are  equal  to 
the  relative  terminal  fall  velocities  of  the  drops,  it  is  possible  to 
express  Dg  in  terms  of  R^  and  Rg  or  p and  Rp  or  even  p and  Rg.  The  intercept 
of  the  function  Dg  with  the  plane  Dg  = O.ly  divides  the  plane  into  regions 
of  coalescence  and  bouncing. 

Fig.  5-4  and  5-5  describe  these  curves  in  the  plane  (p,R^)  for  the 
deformations  Rp  = aR^  and  Rp  = 6Rg/D,  respectively.  Each  deformation 
coefficient  produces  a curve  which  divides  the  plane  into  a coalescence 
region,  below  the  curve,  and  bouncing  above  it. 

If  one  looks  at  a region  of  coalescence  and  reduces  the  size  of  the 
large  drop,  without  changing  p,  there  will  come  a point  at  which  bouncing 
will  start.  This  phenomenon  was  observed  by  many  investigators.  For 
example  Telford  and  Thondike  (1961)  reported  that  drops  smaller  than  35y 
diameter  did  not  coalesce.  Levin  et  al  (1973)  reported  that  bouncing 
occurred  for  all  drops  smaller  than  30-40y  radius.  We  see  that  the  present 
work  predicts  a critical  radius  which  depends  on  the  p ratio. 


Fig.  5 1.  A schematic  graph  representing  the  dependence  of  Dg  on  the  radius 
of  the  small  drop,  Rg,  and  on  the  initial  velocity. 
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Fig.  5-4.  Regions  of  bouncing  and  coalescence  in  the  (p,R^)  plant  for  R^  = aR_. 
The  enclosed  area  is  enlarged  in  Fig.  5-6. 


AREA  OF  COALESCENCE 


The  same  as  in  Fig.  5-A  but  for  R, 
enlarged  in  Fig.  5-7. 


The  enclosed  area  is 


The  results  reported  by  Neiburger  et  al  (1972)  and  Levin  et  al  (1973) 
were  carried  out  in  a situation  in  which  the  two  interacting  drops  moved 
with  respect  to  their  environment.  Hence,  aerodynamic  factors  might  have 
influenced  the  coalescence.  However,  it  is  interesting  to  compare  the 
curves  (p,Rl)  to  their  results.  Figs.  5-6  and  5-7  show  an  enlarged  area 
from  Figs.  5-4  and  5-5.  Neiburger  et  al  (1972)  and  Levin  et  al  (1973) 
observed  no  coalescence  for  droplets  ^ 40p. 

In  Figs.  5-6  and  5-7  x represent  regions  in  which  the  above  authors 
observed  bouncing  and  lines  where  they  observed  coalescence.  It  can  be 
observed  that  the  curve  a = 0.38  in  Fig.  5-6  and  3 = 0.44  x 10®  in  Fig.  5-7 
fit  the  data  well.  These  values  are  very  close  to  the  values  used  in 
the  present  work. 


5-3.  Regions  of  bouncing  and  coalescence  in  non-head-on  collisions 

As  in  the  plane  (p,Rj^)  so  in  the  plane  (0^,V.)  one  gets  curves 
which  separate  between  regions  of  bouncing  and  coalescence.  These  curves 
give  critical  values  of  0^  and  V^. 

i.  Critical  angle  for  coalescence,  C.A.C. 

The  maximum  angle  for  coalescence  is  the  angle  above  which  no 
coalescence  occurs.  In  a similar  way  we  can  define  the  minimum  bouncing 
angle  which  is  the  angle  below  which  no  bouncing  is  observed.  Some  experi- 
ments suggest  that  both  bouncing  and  coalescence  occurs  at  all  angles. 

For  these,  the  above  definitions  are  not, > relevant.  However,  other  experi- 
ments suggested  definite  regions  of  bouncing  and  coalescence.  For  these 
cases  it  is  worth  defining  a critical  angle  of  impact  for  coalescence  - 
C.A.C.  - below  which  coalescence  always  occurs  and  above  which  bouncing 
takes  place.  In  many  experiments  identical  collisions  sometimes  result  in 
bouncing  and  sometimes  in  coalescence.  This  stems  from  the  fact  that  in 
addition  to  the  initial  conditions  other  factors  such  as  waves  on  the 
surfaces,  aerosols  in  the  drops,  and  the  drops'  oscillations  also  influence 
the  coalescence  process.  The  question  then  is  not  what  9.  is  but  rather 
what  is  the  probability  for  coalescence  at  6^.  Under  such  conditions  it 
is  suggested  to  define  an  inpact  angle  for  coalescence,  6,  which  can  be 
expressed: 


TT/2 

sin^  0 = / f(0  ) sin  (20.)d0,  (5-2) 

O 1 11 

where  f(0^)  is  the  probability  for  coalescence  (even  partial  coalescence) 
at  impact  angle  9^. 
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Fig.  5-6.  Regions  of  coalescence  and  bouncing  in  the  (p.^^)  plane  for  ■ aR^.  I 

X and  . represent  bouncing  and  coalescence,  respectively  in  the  j 

data  of  Levin  et  al  (1973).  i 
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sin^0  is  the  collection  efficiency  or  the  coalescence  efficiency  in 
cases  of  impact  on  a stationary  target.  If  a critical  impact  angle  for 
coalescence  exists  than 


and 


f(9^) 


^ 6 

c 


0 > 6 
i c 


sin^0  = /*^sin  6.  d0. 

o i i 


It  seems  that  a critical  inpact  angle  for  coalescence,  CAC,  actually 
exist,  but  that  disturbances  smear  it  out.  The  smaller  these  disturbances 
are  the  closer  the  impact  angle  for  coalescence  comes  to  CAC. 

ii.  Available  data  on  the  critical  Impact  angles 
a.  Impact  on  a stationary  target 

Jayaratne  and  Mason  (1964)  investigated  inpacts  of  drops  on 
horizontal  surfaces  and  found  critical  angles  below  which  coalescence 
occurred  and  above  which  bouncing  was  observed.  They  also  investigated 
inpacts  of  drops  on  curved  surfaces.  For  drops  larger  than  120y  they 
observed,  as  before,  regions  of  coalescence  and  bouncing.  However,  for 
small  drops  they  found  four  regions  of  impact  angles  in  which  coalescence 
and  bouncing  occurred  starting  with  coalescence  at  head-on  collisions. 

Whelpdale  and  List  (1971)  presented  their  results  in  the  form  of 
histograms  from  which  it  can  be  seen  that  coalescence  occurs  at  low  impact 
angles  and  bouncing  occurs  at  higher  angles.  They  estimated  that: 


0 = arc  sin  t — ; — 

c 1 + P 


(5-3) 


b.  Impact  of  two  moving  drops 

Park  (1970),  Nelson  and  Gokhale  (1973)  and  others  have 
observed  that  coalescence  usually  occurs  at  head-on  collisions  and  bouncing 
at  large  inpact  angles.  Many  have  observed  that  in  the  intermediate  angles 
partial  coalescence  usually  occurred. 

Park  (1970)  also  observed  that  the  critical  impact  angles  decrease 
as  the  Impact  velocity  decreases.  A number  of  different  values  of  CAC 
appear  in  the  literature  and  they  vary  from  experiment  to  experiment  (e.g. 

30®  according  to  Ryley  and  Woods,  1965,  1966  and  60®  according  to  Tubman,  1968), 
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Hi.  Critical  in^jact  angle  for  coalescence , CAC,  from  the  present  model 
In  this  section  the  CAC  for  impact  on  a stationary  target  is 

computed. 

In  the  model  as  9.  increases  so  do  x,  (Eq.  4-3)  (Eq.  4-11)  and  D^, 

If  there  is  a 6^  for  which  Dg  •=  O.ly  than  for  lower  0^^  Dg  > O.ly  and 
coalescence  will  occur.  For  larger  6^  Dg  > O.ly  and  bouncing  will  occur. 
This  6^  is  the  CAC.  A computer  program  was  developed  which  provided  the 
CAC  for  given  initial  conditions  Rg»  and  f|), 

Jayaratne  and  Mason  (1964)  and  Park  (1970)  reported  that  CAC  is  a 
monostonic  increasing  function  of  the  impact  speed.  Whelpdale  and  List 
(1971)  found  similar  dependence  for  high  impact  speeds.  However,  for 
lower  Impact  speeds  it  is  difficult  to  conclude  about  the  dependence  of 
®i  '^i  their  resuLts  since  both  increasing  and  decreasing  trends 

are  observed  (see  Fig.  5-9) . In  the  present  model  we  obtain  a monotonlcally 
increasing  functions  of  9j  with  V.  as  is  shown  in  Fig.  5-8  for  Rl  • 250y, 

Rg  * 50y  and  Rg  = 25y.  This  figure  was  drawn  for  Rq  = 0.5  Rg.  For 
Rq  = 0.5  X 10“5  Rg/D  a similar  graph  is  obtained  but  with  different  values. 
According  to  Whelpdale  and  List  (1971)  CAC  is  a monotonlcally  decreasing 
function  of  p (Eq.  5-3).  This  tendency  is  also  obtained  here  (see  Fig.  5-8). 

Whelpdale  and  List  (1971)  presented  their  results  in  histograms  which 
display  the  frequency  of  coalescence  for  each  Aotj^  = 10°,  where  is  defined 
in  Chapter  4 section  a.  Their  results  relate  to  three  values  of  Rj^  and 
four  values  of  (see  Fig.  5-9). 

For  each  there  is  a CAC  within  a range  of  angles  which  we  can  call 
the  critical  range.  The  smallest  angle  in  this  range  (marked  by  a dashed 
line)  is  the  angle  below  which  only  coalescence  occurs.  The  largest  angle 
in  this  range  (marked  by  a dotted  line)  is  the  one  above  which  only  bouncing 
takes  place.  Since  Whelpdale  and  List  (1971)  did  not  provide  a velocity 
spectrum  for  each  case  it  is  difficult  to  con5)are  the  results  of  the  present 
model  to  their  results.  However,  the  use  of  the  two  deformations  (Rp=Q.5Rg 
and  Rp  » 0.5  x 10“®Rg/D)  results  in  CAC  which  falls  within  the  critical 
range  shown  in  Fig.  5-9,  In  fact  the  results  of  the  present  model  fit 
better  as  the  impact  speed  becomes  higher. 

iv.  Critical  impact  velocity 
Available  data 

Tverskaya  (1954)  observed  two  critical  Impact  velocities.  At 
lower  velocities  all  interactions  resulted  in  coalescence  while  at  higher 
velocities  only  partial  coalescence  was  observed.  Impacts  with  velocities 
of  values  in  the  intermediate  range  resulted  in  both  coalescence  and  partial 
coalescence. 


Ilegions  of  bouncing  and  coalescence  in  the  (0.,V  ) plane  for  two 
different  p ratios.  ^ ^ 
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(’omparlson  of  the  present  results  with  the  data  of  Whelpdalc  and 
List  (1971).  Dotted  areas  represent  regions  of  bouncing.  Dashed 
areas  are  regions  of  partial  coalescence.  Regions  in  white  are 
regions  of  coalescence.  Based  on  the  present  work  bouncing  will 
occvir  to  the  right  of  the  dotted  line  and  coalescence  to  the  left 
of  the  dashed  line. 


i.  Definitions 


Coalescence  distance,  x^,  - the  maximum  impact  distance  in 
which  coalescence  occurs. 


Prokhorov  (1949),  (1954)  foimd  that  all  interactions  resulted  in 
coalescence  except  the  ones  at  very  low  Impact  speeds.  Similarly,  Canlas 
(1960)  found  that  bouncing  occurred  at  low  impact  speeds.  As  the  impact 
velocity  increased  and  passed  a critical  value  all  interactions  resulted 
in  partial  coalescence.  When  the  impact  speed  increased  even  further  a 
second  critical  velocity  was  reached  beyond  which  interactions  resulted 
in  coalescence.  Similar  results  were  reported  by  Schotland  (1960). 


Sartor  and  Abbott  (1972)  in  an  experiment  with  two  falling  drops 
observed  a critical  velocity  below  which  bouncing  occurred  and  above 
which  coalescence  took  place.  They  also  observed  a lower  critical  velocity 
below  which  coalescence  occurred. 


As  was  pointed  out  before,  the  curve  in  Fig.  5-8  divides  the  plane 
into  regions  of  coalescence  and  bouncing.  For  each  the  point  on  the 
curve  gives  9^  = 6^,  the  critical  angle.  Similarly,  each  9^  = 9^  on  the 
curve  gives  = V , critical  velocity.  If  is  kept  constant  and  Rg 
is  changed  we  obtain  a curve  (in  the  (9j^,Vj^)  plane)  which  cuts  the 
axis  at  different  points,  this  is  the  critical  velocity  at  head-on-collision. 
We  see  that  the  model  gives  the  critical  velocities  for  each  impact  angle. 


If  a CAC  exists  than  x can  be  written 

c 

X = (R,  + R ) sin  9 (5-4) 

C L S C 

collision  distance,  y^,  - the  maximum  Impact  distance  in  which  collision 
occurs. 


The  relation  between  9 and  V in  the  (G. ,V. ) plane 

c c i’  i 


5-4.  Coalescence  efficiency 


In  this  section  a method  is  suggested  by  which  the  coalescence 
efficiency  and  collection  efficiency  can  be  computed. 


The  relation  between  the  collection  efficiency  E and  x^  is: 


(5-5) 
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The  connection  between  the  collision  efficiency,  and  is 


+ "s)‘ 


(5-6) 


The  relation  between  the  coalescence  efficiency,  E2,  and  is: 


2 


(5-7) 


It  is  easy  to  see  that  the  relation  between  these  equations  is; 


^ = ^1^2 


If  there  is  no  definite  CAC  we  can  define 


= (2/ 


Rs  + h 


f(xi)xidxi) 


where  f(  is  the  probability  that  an  impact  from  an  impact  distance 
results  in  coalescence.  It  follows  that 


and 


X = (R^  + R^)  sine 


li.  Available  data 

The  available  data  from  the  literature  points  out  to  a variety  of 
values  of  coalescence  efficiencies  ranging  from  zero  to  unity.  However  many 
of  the  experiments  expressed  their  results  in  terms  of  full  coalescence 
versus  partial  coalescence  and  bouncing.  Even  the  experimental  part  reported 
on  at  the  beginning  of  this  report  separates  the  two  events.  However,  in 
the  present  theoretical  model  the  coalescence  is  defined  as  any  bridging  of 
the  surfaces  of  drops  regardless  of  whether  it  results  in  full  or  partial 
coalescence. 


>»**:*«< 
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iii.  Coalescence  efficiency  for  impacts  on  a stationary  target 

We  will  define  E„  as  the  coalescence  efficiency  of  drop  impacting 
on  stationary  targets.  In  this  case  E = 1 (see  the.  experiments  of  Whelpdale 
and  List,  1971,  Levin  and  Machnes,  1976,  and  Chapter  2 of  this  report). 

Hence  we  get : 


2R 


= E = 


(■>1  + %y- 


(5-8) 


The  computer  program  which  provides  6 and  x also  computes  E__. 

C C Z K 

We  are  interested  in  values  of  E„  for  inpacts  at  the  relative  terminal 
fall  velocities  of  the  drops.  Tables  h-1  and  5-2  present  E„_  as  a function 
of  and  p for  the  deformations  Rp  = 0.5  R^  and  Rp  = 0.5  x I0“®Rg/D, 
respectively. 

Figs.  5-10  and  5-11  presents  E^^^  as  a function  of  and  p for  the 
above  deformations.  For  clarity,  part  (a) of  the  figures  was  plotted  for 
small  Rj^  and  part(b)for  large  values  of  Rl.  Fig.  5-12  shows  the  relation 
of  E„j-  with  p as  was  suggested  by  Whelpdale  and  List  (1971)  : 


^2R  " (1  + p)2  (5-9) 


The  figure  also  presents  the  results  from  the  present  model  with  both  deforma- 
tions for  the  drop  size  (Rg  = 35y)  used  by  Whelpdale  and  List. 

As  was  already  pointed  out,  Eq.  (5-9)  has  two  main  draw  backs:  i.  There 

is  no  dependence  of  E£j^  on  R^^  (see  Levin  and  Machnes,  1976,  and  Chapter  2 
of  this  report),  ii.  The  lowest  value  of  E2pj  according  to  this  equation 
is  0.25  while  other  experiments  (Levin  et  al,  1973)  found  regions  of  E2j^=0. 
Figs.  5-10  and  5-11  show  that  the  present  model  does  not  have  these  limita- 
tions. In  these  figures  £2^^  is  a monotonir  decreasing  function  of  p to 
the  point  at  which  £2^  = 0.  E.,^  remains  zero  for  larger  values  of  p (it 

should  be  remembered  ^at  the  assumption  in  the  niC'iel  makes  it  most  appli- 
cable for  small  values  of  p) . Both  Whelpdale  and  List  (i971),  and  the 
experimental  results  reported  in  this  report  suggest  similar  dependence 
of  E2j^  on  p. 

Fig.  5-13  describes  E2^  as  a function  of  Rj^  and  p in  a three  dimensional 
graph.  Here  we  see  that  E2j^  increases  as  R^^  increases  up  to  about  lOOy  and 
then  decreases  slowly.  The  intercept  of  the  function  £2^  ” with 

the  plane  (p,Rj^)  gives  the  same  curves  as  in  Fig.  5-3  and  5-4  which  were 
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Fig.  5-10.  The  coalescence  efficiency  as  a function  of  the  ratio  of  the  radii 
P»  for  small  drops  and  (b)  for  larger  drops. 
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.escence  efficiency  as  a function  of  the  p ratio  for 
small  drops  and  (b)  for  larger  drops. 


obtained  from  head-on  collisions. 


If  Fig.  5-13  and  2-11  are  compared  one  sees  similar  trends  in  the 
range  R,  > 125y  in  which  the  experiments  were  carried  out.  Looking  at 
Fig.  2-10  one  can  see  that  experimentally  some  increase  in  £2^^  with  an 
increase  in  was  observed  in  qualitative  agreement  with  the  left  side 
of  Fig.  5-13  and  part  a of  Figs.  5-10  and  5-11. 

In  Tables  5-1  and  5-2  and  Figs.  5-10,  5-11  and  5-13  values  of  £2^ 
were  calculated  for  Vj^,  the  relative  terminal  velocities  of  the  drops. 

The  model  is  also  capable  of  calculating  these  efficiencies  for  any 
impact  velocity. 

iv.  Coalescence  efficiency  of  moving  drops 

Behind  the  definitions  of  collision  and  coalescence  efficiencies 
rests  the  assumption  that  the  approach  of  two  drops  can  be  divided  into 
a collision  region  in  which  aerodynamic  effects  dominate  and  a coalescence 
region  in  which  coalescence  effects  dominate. 

If  p is  small  than  when  the  deformation  begins,  the  small  drop  will 
already  be  under  the  influence  of  the  boundary  layer  of  the  larger  drop 
(see  Whelpdale  and  List,  1971)  and  the  present  model  is  applicable. 
However,  if  p is  larger  the  coalescence  effects  would  affect  the  closer 
surfaces  of  the  drops  while  the  further  surfaces  would  still  be  affected 
by  the  free  stream  flow.  For  such  a case  the  separation  between  collision 
and  coalescence  is  not  fully  justifiable. 

Let  r^  be  the  point  at  which  the  coalescence  effects  start  to 
dominate.  Then  if  a drop  starts  its  collision  at  r > rj,  and  initial 
conditions  (Rl,  Rs,  V^,  x.,n)  it  will  arrive  at  Vq  with  conditions 
(Rt  , Rg  , Vq,  3fg,r|)  as  a result  of  the  collision  effects.  At  r^  4 Vq 
ana  ¥ Xq-  The  present  model  follows  the  approach  from  r^  to  Dg, 
as  was  discussed  in  the  last  three  chapters.  Therefore,  values  of 
coalescence  efficiencies  for  these  cases  can  be  computed. 


CHAPT£R  6.  CONCLUSIONS 


This  report  summarises  results  of  the  coalescence  efficiency  of 
colliding  drops  which  were  obtained  from  laboratory  experiments  and  a 
theoretical  model.  Both  theory  and  experiment  suggest  that  the  coale- 
scence efficiency  is  a function  of  both  drop  radii  and  drop  size  ratio. 

£ven  though  no  clear  dependence  of  coalescence  efficiency  on  impact 
speed  was  observed  experimentally,  the  theory  does  predict  some  depen- 
dence on  velocity.  Generally  it  is  observed  that  the  coalescence  efficiency 
decreases  with  increase  in  the  size  of  the  large  drop,  Rj^,  and  the  drop 


ratio,  p.  However,  for  small  a decrease  in  the  efficiency  was  also 
shown. 

Both  theory  and  experiment  demonstrate  the  strong  dependence  of 
coalescence  on  iin)act  angle.  The  theory  is  able  to  predict  the  efficiency 
of  coalescence  for  different  impact  speeds  and  not  only  the  relative 
terminal  velocities  of  the  drops.  It  is  also  capable  of  calculating  the 
coalescence  efficiency  of  all  Impacts  on  stationary  targets  and  the 
coalescence  efficiencies  for  moving  drops  when  the  p ratio  is  small. 

From  the  experiments  two  impirical  equations  were  obtained  which  relate 
the  coalescence  efficiencies  and  partial  coalescence  efficiencies  with 
p and  IL  . These  could  easily  be  used  in  numerical  models  of  cloud  and 
rain  dewlopment.  Even  though  the  present  results  provide  physical 
understanding  and  numerical  values  of  the  coalescence  efficiency,  the 
problem  is  far  from  being  fully  understood.  It  is  hoped  that  this 
investigation  will  continue  in  order  to  find  values  of  coalescence 
efficiency  for  much  smaller  drops  and  for  all  p ratios. 

It  should  be  pointed  out  that  the  theoretical  derived  values  of 
the  coalescence  efficiency  were  generally  lower  than  our  experimental 
results. 


We  feel  that  the  inclusion  of  surface  waves  on  the  deformed  surfaces 
in  the  theory  would  Increase  the  theoretical  values  and  bring  them  closer 
to  those  of  the  experiment.  We  plan  to  continue  our  theoretical  investi- 
gation in  this  direction. 
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